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ABSTRACT 

The optical and X-ray light-curves of long Gamma Ray Bursts (GRBs) often show a com- 
plex evolution and in most cases do not track each other This behaviour can not be easily 
explained by the simplest standard afterglow models. A possible interpretation is to consider 
the observed optical and X-ray light-curves as the sum of two separate components. This sce- 
nario requires the presence of a spectral break between these bands. One of the aims of this 
work is to test whether such a break is present within the observed Swift XRT energy range. 
We analyse the X-ray afterglow spectra of a sample of 33 long GRBs with known redshift, 
good optical photometry and published estimate of the host galaxy dust absorption Ay"**. We 
find that indeed in 7 bright events a broken power-law provides a fit to the data that is better 
than a single power-law model. For 8 events, instead, the X-ray spectrum is better fitted by 
a single power-law. We discuss the role of these breaks in connection to the relation between 
the host hydrogen column density N^°^^ and ^y°''' and check the consistency of the X-ray 
spectral breaks with the optical bands photometry. We analyse the optical to X-ray spectral 
energy distributions at different times and find again consistency with two components inter- 
pretation. 
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1 INTRODUCTION 

The fast re-pointing capabilities of the Swift satellite (Gehrels et 
al. 2004) allowed to reveal the early time afterglow behaviour of 
Gamma-Ray Bursts (GRBs) and its unforeseen complexity. Sev- 
eral interpretation have been proposed to account for the unex- 
pected light-curves evolution. This is often characterised by an 
early time steep flux decay starting after the end of the gamma-ray 
prompt emission, followed by a shallower and a subsequent steeper 
decay, the latter corresponding to the typical afterglow observed in 
the pieSwift era (see e.g. Nousek et al. 2006, Zhang et al. 2006). 

Great efforts have been made to explain the origin of the shal- 
low decay phase which can last from some hundreds up to hundred 
of thousand seconds as it cannot be explained in the frame of the 
simplest "standard" forward shock fireball models. In Ghisellini et 
al. (2009) (hereafter G09), we presented a summary of the proposed 
interpretations (see also Zhang 2007 for a review). 

In the last years, the increasing number of well sampled light- 
curves allowed to examine simultaneous optical and X-rays after- 
glow light-curves of several long GRBs (e.g. Curran et al. 2009, 
De Pasquale et al. 2009). In G09 we analysed the broad band op- 
tical to X-ray rest frame temporal behaviour of a sample of 33 
GRBs observed by Swift XRT with known redshift, published host 
galaxy dust absorption estimate and good quality optical follow up. 
In some cases the optical and X-ray temporal evolution are very 
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different. We proposed that the light curves behaviour is due to 
the sum of two separate components. The first one is assumed to 
originate from a standard forward (external) shock, as described 
by Panaitescu & Kumar (2000). The second component is treated 
in a purely phenomenological way with the aim of minimising the 
number of free parameters. A possible physical origin for it can be 
provided within the so called "late prompt" scenario described by 
Ghisellini et al. (2007). In G09 we found that this two component 
modelling is able to well reproduce all of the optical and X-ray 
light curves of the GRBs of the sample (once the early steep decay 
phase and the flaring activity that sometimes appears in the X-ray 
light-curves are excluded). 

In order to test the consistency of the two components inter- 
pretation and make a first step towards a more physical scenario, 
it is important to verify whether the observed X-ray spectra and 
the optical to X-rays spectral energy distributions (SEDs) are in 
agreement with what predicted by the light-curves modelling. 

In a scenario where the optical and the X-ray emission are 
due to different processes, the component accounting for the opti- 
cal spectrum has to break in order not to dominate also in the X-ray 
band. Conversely a break to a harder spectral index is also required 
towards the soft end of the X-ray component, not to interfere with 
the optical emission. However, a spectral break (e.g. the cooling 
break frequency of the synchrotron emission mechanism) between 
the optical and X-ray bands is sometimes expected also in the stan- 
dard afterglow scenario (see e.g. Figs. 10 and II in Nardini et al. 
2006). In both scenarios the spectral break can be at frequencies 
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within tlie observed XRT 0.3-10 keV energy range. If this is the 
case, spectral fits provide not only the break frequency but also the 
spectral slope below it. This additional information makes these 
events the best candidates to test the two component light-curve 
modelling from a spectral point of view. If the optical and X-ray 
light-curves are dominated by the same component, the observed 
optical fluxes must be consistent with the extrapolation of the X- 
ray low energy spectrum. If the light-curves are instead dominated 
by different components, the X-ray spectrum extrapolation should 
not significantly contribute to the observed optical flux (see ^for 
a more detailed discussion). 

In this work we analyse the XRT spectra of the GRBs in the 
G09 sample to check for the presence of such a break. In order to 
test whether the X~ray spectral break is consistent with what seen 
(simultaneously) in the optical, in 21 we examine the optical to X- 
ray SEDs sampled at different times along the light-curves. Such 
a combined analysis of broad band light-curves and optical to X- 
ray SEDs represents a crucial consistency check for our proposed 
interpretation. 

As discussed in section i]3.3l an interesting outcome of the 
spectral analysis is related to the apparent "discrepancy" between 
the amount of the X-ray absorption (as measured by the hydrogen 
column density A^'h"'"') and the optical extinction A^''^ in the host 
frame. The value of A^^°''' inferred from fitting the X-ray spec- 
tra with a single power-law is often at odds (for standard gas-to- 
dust conversions) with the relatively small Av""* evaluated through 
the analysis of the optical SEDs (see e.g. Galama & Wijers 2001; 
Stratta et al. 2004; Kann, Klose & Zeh 2006, Schady et al. 2007). 
If the intrinsic X-ray spectrum can be well modelled by a broken 
power-law, then the required N^'^^ is smaller than what required 
by a single power-law fit), ameliorating the 7Vh°^'~^v°^' disagree- 
ment. 



2 THE SAMPLE 

The sample comprises the 33 long GRBs considered G09, whose 
selection criteria were: the knowledge of the GRB redshift, a good 
photometric coverage. Swift XRT observations and a published es- 
timate of the host galaxy dust absorption A^^^ . When different 
values of Ay""* are reported in the literature we chose the estimate 
derived from a direct analysis of the optical spectral energy dis- 
tribution rather then that obtained by a combined analysis of the 
optical to X-ray SEDs. If only the latter is available we discuss the 
effects of possible alternative solutions, through a direct analysis of 
the SEDs (see ©. 



3 XRT DATA REDUCTION AND SPECTRAL ANALYSIS 

We analysed the XRT data of the events in the sample with the Swift 
software package v2. 9 distributed with HEASOFT {v6.6). The XRT 
data were reprocessed with the XRTPIPELINE toof^. The spectra 
were extracted in both WT and PC mode with the standard grade, 
applying, when required, the correction for pileup (Moretti et al. 
2005, Romano et al. 2006; Vaughan et al. 2006). The extraction 
was in boxes (WT mode) or circular regions (PC mode) of typi- 
cal widths as discussed in Evans et al. (2009). Background spec- 
tra were extracted in same-sized regions far from the source. For 

1 Part of the XRT software, distributed with HEASOFT package: 
|http : / /heasarc ■ gsf c ■ nasa ■ gov/heasof t/| 



all of the spectra we created Ancillary Response Files with the 
xrtmkarf tool and used the calibration database updated to De- 
cember 2008. The spectra were re-binned in order to have a min- 
imum of 20 counts per energy bin (15 for the faintest events) and 
energy channels below 0.3 keV and above 10 keV were excluded 
from the analysis. The XSPEC(vii.5.2) software was utilised for 
the analysis. For bursts with particularly bright X-ray emission we 
also performed a time resolved spectral analysis in order to check 
for the possible spectral evolution. Since we are not considering 
XRT data that are simultaneous to the BAT 7-ray detection, the 
steep early time phase and the flaring activity are not considered. 

3.1 Single absorbed power-law model 

Following the conventional analysis of X-ray GRB spectra we fit- 
ted all the spectra with a model composed by a power-law with two 
absorption components at low energies, wabs and zwabs. The 
first one corresponds to Galactic absorption and its column den- 
sity A'^^''' is fixed to the Galactic value (from Kalberla et al. 2005). 
The second absorption is due to the material located at the redshift 
of the source and its column density A'^^°'^' is let free to vary. The 
90% confidence intervals on the best fit parameters are obtained 
with the error command in XSPEC. All the spectra returned a 
good fit with such a model, with x^/dof close to unity. The best 
fit parameters are in a good agreement both with the results of the 
automatic XRT data analysis tool available on lintQ developed by 
Evans et al. (2008, 2009) and with the values reported in the litera- 
ture (summarised also in Tab. 1 of G09). The results of the fits are 
reported in Tab. |3.1| 

3.2 Broken power-law model 

In order to test for the presence of possible spectral breaks within 
the XRT energy range we selected the GRBs whose spectra have 
high signal-to-noise, namely those which, after the applied rebin- 
ning, had a minimum of 50 energy bins. This choice, on average, 
corresponds to a minimum of 1000 counts per spectrum. We found 
20 events fulfilling this condition. In the excluded 13 cases (i.e. 
GRB 050319, GRB 050408, GRB 050525A, GRB 050801, GRB 
050824, GRB 05 1 1 1 1, GRB 060512, GRB 060526, GRB 060904B, 
GRB 060927, GRB 070125, GRB 071010A and GRB 080310) the 
spectrum in the considered time intervals has too low S/N for fitting 
a broken power-law model which has two more free parameters (i.e 
the spectral index of the second power-law component and the en- 
ergy break between the two power-laws) with respect to the single 
power-law model with galactic and intrinsic absorption. 

We used 2 absorption components for the broken power-law 
models, as described in the Section [3711 The break energy Eh be- 
tween the low and high energy power-laws spectral indices (/3x,i 
and /3x,2, respectively) was left free to vary in the 0.3-10 keV 
energy range. Clearly a significant broken power-law fit should 
result in statistically different /3x,i and /3x,2. Therefore no pre- 
determined relation between the model parameters was assumed 
(as done for instance if the emission process is assumed to be syn- 
chrotron (/3x,i = /3x,2 — 0.5). 

The broken power-law with a free rest frame A^^"*^' model 
(hereafter ABP) has 5 free parameters while the absorbed single 
power-law model (hereafter AP) has 3 free parameters that are a 

^ |http : / /www ■ swift ■ ac ■ uk/ xrt-spectra/| 
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Table 1. Results of the single power-law fitting. For each GRB we report: the redshift, the time interval in which the spectrum was extracted, the unabsorbed 
spectral index /Jx, the hydrogen column density at the host Af^°^' , the reduced and number of degrees of freedom, the host galaxy visual extinction 
taken from the literature, and the references for redshift and A^"^*. References: Ber05a: Berger et al. (2005a); Sti05: Still et al. (2005); Fyn05a: Fynbo et al. 
(2005a); Kan08: Kann et al. (2009); Fyn06: Fynbo et al. (2005b); Wat06: Watson et al. (2006a); Ber05b: Berger et al. (2005b); dUP07: de Ugarte Postigo 
(2007); Cen05: Cenko et al. (2005); Hol07: Holland et al. (2007); Fol05: Foley et al. (2005); Che05: Chen et al. (2005); Sta05: Starling et al. (2005); DeP07: 
de Pasquale et al. (2007); Fyn05d: Fynbo et al. (2005c); Sch07: Schady et al. (2007); Pro05: Prochaska et al. (2005); Fyn05d: Fynbo et al. (2005f); Jak06: 
Jakobsson et al. (2006); Hil05: HUl et al. (2005); Cen06b: Cenko et al. (2006b); Mis07: Misra et al. (2007); Fyn06: Fynbo et al. (2006a); Cuc06: Cucchiara 
Fox & Berger (2006); Cur07: Curran et al. (2007); Pro06: Prochaska et al. (2006); E1106: Ellison et al. (2006); BI0O6: Bloom et al. (2006); Ber06: Berger & 
Gladders (2006); Tho08: Thone et al. (2008) Pri06: Price Berger & Fox (2006); Man07: Mangano et al., (2007); Th606: Thone et al., (2006); Gru07: Grupe 
et al. (2007); Fug06: Fugazza et al. (2006); R0IO6: Rol et al. (2006); Fyn06b: Fynbo et al. (2006b); RuV07: Ruiz-Velasco et al., (2007); Osi06: Osip Chen 
& Prochaska (2006); Mun07: Mundell et al. (2007); BI0O6: Bloom Perley & Chen (2006), Pag07: Page et al. (2007); Per()8a: Perley et al. (2()08a); Jau07: 
Jaunsen et al. (2007); Tro()7: Troja et al. (2007); Fox et al. (2007); Per()7: Perley et al. (2007); Per08b: Perley et al. (2()08b); Pro07: Prochaska et al. (2007); 
Cov08: Covino et al. (2008); Pro08: Prochaska et al. (2008); PeB08: Perley & Bloom (2008a); Vre08: Vreeswijk et al. (2008); Blo09: Bloom et al. (2009). 



subset of the ABP model ones, models are nested with a progres- 
sion of 2 free parameters so an ABP model fitting is considered an 
improvement of the AP model one if Ax^ = Xap — Xabp > 4.6 
(90% confidence). A similar choice was also done by Butler & 
Kocevski (2007): they considered as acceptable a more complex 
model (with an additional free parameter) if Ax^ > 2.7. 

In 7 events (i.e. GRB 050802, GRB 050820A, GRB 060210, 
GRB 060729, GRB 061007, GRB 061126, GRB 080319B) the fit 
with the ABP model restilted in an acceptable x^/dof and the 5 
free parameters of the ABP model were constrained with accept- 
able uncertainties (i.e. a minimum is found inside the parameters 
definition range also considering their uncertainties). Usually both 
the high energy photon index (/3x,2) and E\, are well constrained 



( typical errors of about 0.1 for the spectral index /3x,2 and 0.15 
keV for the break energy £b while /3x,i) and iV^""' are affected by 
larger - but still acceptable uncertainties (about 0.2 and 50%, re- 
spectively) (see Tab. 2). For all the 7 events the improvement of the 
ABP fit with respect to the AP one yields at least a 90% significant 
improvement. 

In 8 cases (i.e. GRB 050318, GRB 050401, GRB 050416A, 

GRB 050922C, GRB 060614, GRB 060908, GRB 070110, GRB 
071003) the ABP model is not preferred to the AP one, either be- 
cause /3x,i is equal to /3x,2 within their errors or E\, results outside 
the considered energy range. 

In 5 GRBs (i.e. GRB 060124, GRB 060206, GRB 061121, 
GRB 050730, GRB 060418) although the ^ of the ABP model is 



© 0000 RAS, MNRAS 000, 000-000 



4 M. Nardini, G. Ghisellini, G. Ghirlanda, A. Celotti 



lower than that of the AP model, the improvement of the fit is not 
statistically significant. 

We re-analysed the spectra of these 5 events assuming an ABP 
model with 4 parameters, namely with N^^^ frozen to the value 
estimated from A^""^ assuming that the A^^^ - iVn"''* relation re- 
ported by Schady et al. (2007) (their Eqs. 1, 2 or 3). For each burst 
we choose the conversion corresponding to the extinction curve 
adopted to obtain the A^^^ from the analysis of its optical SED. 

In Tab. 3 we report the best X-ray spectral fit parameters val- 
ues for these 5 events and the related Xr- For all of these events we 
obtain a good fit to the data with Xr values close to unity like in 
the AP case. 

As the AP parameters are no more a subset of the parame- 
ters of this model (i.e. they are not nested models) the Ax"^ does 
not provide statistical information on the fit improvement (e.g. Pro- 
tassov et al. 2002). 

No a priori relation between /9x,i and I3x,2 was assumed 
and the extremely hard /9x,i obtained for GRB 061121 and GRB 
060418 cannot be easily accounted for by the standard emission 
processes. Given the uncertainties in the inferred A'^^°'^', we then 
fixed or constrained the value of in two ways, assuming: i) 
the relation A/3 = /3x,i — Px,2 ~ 0.5; ii) Px,i = 0. In both cases 
the best fit returns the same x^d value. For GRB 061121 the de- 



rived columns are Ny, 



and iV, 



0.58 



host 



+o:i3 X 10^2 cm-2 (AP = 0.5) 



GRB 060418 A^h""* 



O.Utoil X 10^2 cm"^ (Px.i = 0), while for 
•28to.?6 X 10^2 cm"^ (A/3 = 0.5) and 
jynost Q 25 ^ iQ^2 cm"^ (Px,i = 0). We conclude that for 
these two bursts the data cannot robustly constrain the low energy 
spectral slope as an acceptable fit can be obtained for not so ex- 
treme values of /3x,i- The column densities obtained in these cases 
are intermediate between the ones obtained through the ABP and 
the AP models, in agreement with what found by Schady et al. 
(2007) when assuming A/3 = 0.5. 



3.3 Discussion on the X-ray spectral analysis 

The breaks that we have found are all in the range between 0.6 and 
1.6 keV. This is likely due to the fact that the peak of the effective 
area of the Swift-XRT is 1.5 keV and it is therefore easier to find 
a spectral break when it falls around this energy. A break at < 
0.6 keV is hardly detectable and therefore we cannot exclude its 
presence in most of spectra. Also a break at E^, > 3 keV can not 
be easily detectable with the available late time X-ray spectra but 
we expect in this case to obtain a quite hard spectrum fitting with 
an AP model (i.e. a single power law spectral index with a value 
similar to the obtained /3x,i). Since the values of /3x reported in 
Tab. l3.1l are usually not so hard we do not expect we are missing a 
large number of > 3 keV breaks. 

It has been already pointed out (e.g. Galama & Wijers 2001, 
Stratta et al. 2004, Schady et al. 2007, Starling et al. 2007, Watson 
et al. 2007) that the A^h"''* derived from the fit of the X-ray spectra 
are usually quite large (up to a few xlO^^ cm~^). On the other 
hand the values inferred for host galaxy dust reddening A\^^^ are 
usually small. This is true also for the GRBs in our sample as shown 
in Fig.[T] 

From the observational point of view the large A'^^°'^' derived 
from the fitting of the X-ray spectrum corresponds to a deficit of 
counts below approximately 1 keV with respect to the extrapola- 
tion of a single power-law model. In principle, this deficit could 
instead be due to an intrinsically curved or a broken power-law 
spectrum. For the 7 GRBs for which the ABP model gave a better 
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Figure 1. Rest frame column densities N^°^^ (obtained from fitting a sin- 
gle power-law model to the X-ray data) versus the visual extinction 
in the GRB host galaxy for all 33 GRBs of the sample. The thi'ee curved 
lines correspond to the Afjj.x versus Ay relations observed in the Milky 
Way and in the Small Magellanic Cloud as described by Eqs. 1, 2 or 3 in 
Schady et al. (2007). 



fit (with respect to the AP model) we can check if the obtained val- 
ues of A^g°^* are in agreement with the optical extinction assuming 
a A^H.x vs. Av relation (e.g. Eqs. 1, 2 or 3 in Schady et al. 2007). 

Fig. [2] shows the values of Ar^°=* obtained with the ABP 
model fitting versus (filled circles). For comparison, also the 

N^°''^ values obtained with the AP model fitting (empty squares) 
are reported. The solid lines represent the Milky Way and Small 
Magellanic Cloud like relations as in Fig.[T] For 5 GRBs the un- 
certainties on 7Vh°^* are quite large, making these values consistent 
with zero, i.e. they must be considered as upper limits. These limits, 
always smaller than the N^°^^ values obtained with the AP model, 
are consistent with the observed A\°''^. For the remaining 2 GBRs 
(GRB 060210 and 080319B) however the value of 7V^5°"* are still 
somewhat larger than what expected by the standard gas-to-dust 
relation, though clearly the disagreement with is less pronounced. 

While the presence of an intrinsic break in the emitted X-ray 
spectrum can solve or mitigate the problem of an excess of A'^^"''' 
with respect to the optical reddening for a fraction of events, this 
can not be considered as a general solution of this issue, on the 
basis of different indications. 

As the excess is observed in a large fraction of GRBs, this 
would imply that the observed X-ray spectrum is almost always a 
broken power law, with a break in the rather narrow 0.5-1.5 keV 
energy range, even if the redshifts of these bursts are different. 

Furthermore we can directly exclude the presence of a spec- 
tral break inside the observed XRT spectrum for about half of the 
analysed events. In general, these events have an intermediate/high 
^H°''' (when fitted with the AP model; see Fig.[3]compared to Fig. 

13. 
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GRB 


z 


1021 cm-2 


/3x,i 


keV 


/3x,2 


Xr (dof) 


prob 




050802 
050820A 
060210 
060729 


0.55 
2.612 
3.91 
0.54 


0.6±0.6 

9 7+2.2 

^■^-2.2 

4+7 

+■-4 

0+0.2 


58+!;;3 

~^ — 0.14 
U.OJ_Q 20 

"■-''^-0.22 
(,^0+0.24 
'^■-'-'-0.09 


1 64+0-«3 

— 0.64 
1 05+0-^0 
l-"-'-0.33 
■ ,,+0.23 
'■'■'-0.17 
1 1T+0.13 
'■'■'-0.10 


0.95±0.12 
1.00±0.07 
1.12±0.07 
1.04±0.04 


0.99 (157) 
0.947 (141) 
0.99(183) 
0.97 (288) 


6.5e" 
7.7e- 
8.5e- 
2.9e- 


-3 
-2 
-2 
-3 


061007 


1.26 


^■'-'-0.9 
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3.9e- 


-4 
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1.1588 
0.937 


1 7+3.6 

(,7+0.2 
'-'■'-0.2 


-0.16+?f 
i'-+y-o.io 




0.74±0.08 
0.81±0.01 


1.056(141) 
1.27 (608) 


9.8e- 
8.6e- 


'2 
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Table 2. Results of the fit to the X-ray spectra with the absorbed broken power-law for the 7 bursts for which the model parameters are constrained. Note that 
P represent the energy spectral index (/3 = F — 1). The analysed spectra have been extracted in time intervals as in the third column of Tab. IxT] 



GRB 
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keV 
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X\ (dof) 


050730 
060206 
060124 
060418 
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1.22 (94) 
0.98 (86) 
0.81 (106) 
0.85 (90) 
0.88 (120) 



Table 3. Results of the absorbed broken power-law model fitting obtained by freezing the value of N^'^^ to that estimated from A^™* through Eq. 1, 2 or 3 
in Schady et al. (2007) (see text). The analysed spectra have been extracted in time intervals as in the third column of Tab. IXT] 
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Figure 2. Rest frame column densities Af^™* versus visual extinction 
y^hoat fo,. the 7 GRBs in which a broken power-law model gave an ac- 
ceptable fit (see text). Filled circles represent the column densities obtained 
from an absorbed broken power law fit to the XRT spectra with the local 
absorption fixed to the Galactic N^'^^ values while empty squares indicate 
the Af^°=' obtained from a single power-law fitting for the same events. 
The three curved lines represent the A^h.x versus Ay relations observed in 
the Milky Way and in the Small Magellanic Cloud as described in Schady 
et al. (2007). 



Figure 3. Rest frame column densities Nj^°'^^ versus the visual extinction 
^host obtained from the AP model for the events where the ABP model 
is excluded. The two curved lines represent the A'^h.x versus Ay relations 
observed in the Milky Way and in the Small Magellanic Cloud as described 
in Schady et al. (2007). 

4 OPTICAL TO X-RAYS SEDS 

In this section we present the broad band SEDs for the 7 events in 
which we found a break in the X-ray spectrum to check whether the 
spectra, when extrapolated to lower frequencies, are consistent with 
the available optical photometry at the same epoch. Then the optical 
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to X-ray SED information will be combined with the results on the 
decomposition of the light curves behaviour in two components, as 
suggested by G09. 

A key test of the two component light curve modelling by G09 
is to verify whether the spectral properties at different times are 
consistent with what inferred by the light curve de-convolution. The 
7 GRBs whose X-ray spectrum is consistent with the presence of a 
break in the X-ray band are the best candidates to perform the test 
as shown in fig.|4] 

For each burst we select epochs where simultaneous optical 
photometry and XRT observations are available in order to use the 
most complete spectral information available and to avoid (if possi- 
ble) flux extrapolations. This limits the number of optical to X-ray 
SED considered. 

The X-ray spectrum is extracted from a time interval around 
the selected epoch in order to have at least 50 energy bins and it is 
re-normalised to the 0.3-10 keV flux obtained from the light curve. 
The spectral index plotted in the SEDs are the ones reported in Tab. 
2. As done in G09 we used the light curves from the Swift reposi- 
tory (see Evans et al. 2009). The counts to flux conversion instead 
was done using the values from our broken power law spectral fits. 
When the optical and X-ray bands light curves track each other, 
i.e. they are dominated by the same component, one single SED is 
considered; when instead they show different temporal behaviours 
we considered more SEDs to test the modelling at different phases 
of the evolution of the two components. 

In the following we present the results for 6 of the 7 GRBs sep- 
arately. The complexity of the remaining one (i.e. GRB 0803 19B) 
prompted us to discuss it in the details in a dedicated paper (Nar- 
dini et al. 2009 in preparation). We anticipate that no event shows 
an optical to X-rays SED that is inconsistent with the presence of 
a break in the XRT band and the two component interpretation. 

4.1 GRB 050802 

The optical light curve photometric data are mainly from Gates 
et al. (2007) together with later time R band data from GCNs 
(Pavlenko et al. 2005, Fynbo et al. 2005). The Swift UVOT filters 
uvm2 and uvw2 are strongly affected by Lya dumping and are not 
considered. The X-ray light curve has been modelled as the com- 
bination of "standard afterglow" emission dominating at early and 
at late times (before 700 s and after about 10 ks, rest frame) and 
"late prompt" emission dominating in between. The "standard af- 
terglow" component instead describes the evolution of the optical 
flux during the whole period of the follow up (see Fig.[5](- 

We extracted the optical to X-ray SED around an observed 
time of 2500s (920 s rest frame) when the optical light curve is 
dominated by the standard afterglow and "late prompt" emission is 
becoming predominant in the X-rays. The X-ray spectrum repre- 
sented in Fig. [5] has been extracted in the time interval reported in 
Tab. 1. Schady et al. (2007) estimated a non negligible host galaxy 
dust absorption (yly"'"' = 0.55 ± 0.1) on the basis of a Milky 
Way extinction curve and assuming a power-law spectrum con- 
necting the optical and X-ray bands. By considering the optical 
bands alone we find a similar Av°°' =0.6 with an optical spectral 
index /3o ~ 0.9. 

The SED, plotted in Fig.jS] is consistent with the optical and 
the X-ray emission being dominated by different components (note 
that also Gates et al. 2008 and de Pasquale et al. 2009 found a simi- 
lar inconsistency) with a spectral break falling in the observed XRT 
energy range, as indeed obtained from the X-ray spectral analysis. 
Note that the X-ray spectra shown in this section have been "de- 
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Figure 4. Sketch illustrating the possible different cases for the relation be- 
tween light curves behaviour in terms of the two component decomposition 
and SEDs. The left panels refer to the X-ray (upper curves) and optical 
(lower curves) light curves, while the corresponding expected optical to X- 
ray SED are shown in the right-hand panels. The bottom right panel shows 
the standard "afterglow-afterglow" case, i.e. both light curves are domi- 
nated by the afterglow component, with a cooling break appearing first at 
X-ray energies. The vertical grey line indicates the time of the extraction of 
the SED. v\, and Vc represent the break and cooUng frequency, respectively. 

absorbed" both for the galactic and the host frame (when present) 
contributions. 



4.2 GRB 050820A 

The photometric data are from Cenko et al. (2006a, 2009). 
SwiftrBPS triggered on a precursor about 200 s before the main 
event (Cenko et al 2006a, Burlon et al. 2009). Our reference time 
is set at the trigger time and we do not consider the prompt X-ray 
emission detected before the end of the main 7-ray event. The X- 
ray light curve is dominated by the "late prompt" component up 
to 200 ks (720 ks in the observer frame), and by a "standard af- 
terglow" component after then. The "standard afterglow" emission 
instead prevails during the entire duration of the optical light curve 
but in a time interval around 5 ks (18 ks observer frame) where its 
contribution becomes comparable with the "late prompt" one (see 
Fig.|6j. 

We extracted two SEDs in order to test the modelling at two 
different light curve phases. The first one at about 20 ks in the ob- 
server frame (~ 5500 s rest frame) where the "late prompt" gives 
the maximum contribution in the optical light curve and the avail- 
able photometry is richer (7c, Rc, V, g and B bands). Cenko et 
al. (2006a) estimated a /3c, — 0.77 with negligible host galaxy 
dust absorption while Kann et al. (2009) inferred an y4v°°' = 
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Figure 5. Top panel: X-ray (in grey) and optical (different symbols, as 
labelled) light curves of GRB 050802 (in the rest frame time). Lines indicate 
the model fitting: afterglow component (dashed line), "late prompt" one 
(dotted line) and their sum (solid line). Black lines refer to the X-rays, light 
grey (red in the electronic version) to the optical. The vertical line marks the 
time at which the SED is extracted. Bottom panel: Optical to X-ray uF,y 
SED at about 2500 s (observed time, corresponding to 920 s rest frame) 
after trigger. 

The dashed and soHd Hues show respectively the best fit (with the ABP 
model) to the X-ray spectrum (the spectral parameters are reported in Tab. 
2) and the uncertainties on the slope of the low energy spectral index /3x,i- 
The dotted line shows the best fit (with the AP model) to the X-ray spec- 
trum. 

0.065 ± 0.008. We used the latter estimate and obtained /3o « 0.7. 
This first SED is plotted in Fig. [6] and shows that the optical flux 
lies slightly above the extrapolation of the broken power law that 
best describes the XRT spectrum, but as the uncertainties on /3x,i 
are quite large the optical flux is fully consistent with the extrapola- 
tion. In this SED the X-ray data are extracted from the time interval 
reported in Tab. 1 . 

As mentioned at 5500 s (rest frame) the optical flux is due 
to a similar contribution of the "standard afterglow" and the "late 
prompt" component. The the cooling frequency is already redward 
of the considered optical bands and the "standard afterglow" has 
/3o ~ 0.92 (corresponding to an emitting particle distribution with 
slope p = 1.85). In the "late prompt" component modelling, the 
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Figure 6. Top panel: X-ray and optical light curves of GRB 050820A (rest 
frame time). Same notation as in Fig. |3] The thin vertical line represents 
the jet break time with its estimated errors (see Ghirlanda et al. 2007 and 
references therein). The thick vertical lines mark the times at which the 
SEDs are extracted. Middle and bottom panels: Optical to X-ray uF,^ SED 
around 20 ks (middle) and 150 ks (bottom) after trigger in the observer 
frame (corresponding to 5.5 ks and 41.5 ks in the rest frame). The dashed 
and solid lines show respectively the best fit (with the ABP model) to the 
X-ray spectrum (the spectral parameters are reported in Tab. 2) and the 
uncertainties on the slope of the low energy spectral index /3x i . The dotted 
line shows the best fit (with the AP model) to the X-ray spectrum. 
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low energy spectral index is instead /3o = 0.45 (see Eq. 3 in G09), 
consistent within errors with /3x,i- The intermediate value of the 
observed optical slope is thus consistent with the predictions of the 
two component modelling. 

We considered a second SED at about 150 ks after the trig- 
ger (observer frame, corresponding to 41 ks rest frame). The plot- 
ted X-ray data are from the time integrated spectrum of the com- 
plete second XRT observation. In this phase the X-ray light curve 
is dominated by the "late prompt" while the "standard afterglow" 
dominates the optical emission. The combined SED is plotted in the 
bottom panel of Fig.|6]and confirms the proposed scenario: the op- 
tical data are at this time much brighter than what predicted by the 
extrapolation (with slope /9x,i) of the X~ray spectrum to the optical 
bands . Even though at these late times there are only 3 available 
photometric points (Jc, Rc and V band) and the V band flux is af- 
fected by a large error, the optical SED is well fitted by a softer 
/3o = 0.95 that is closer to the value predicted for the "standard 
afterglow" component. 

4.3 GRB 060210 

The optical afterglow has been observed in the R and / bands (pho- 
tometric data from Curran et al. 2007), while, because of the high 
redshift, smaller wavelengths bands are not observable due to the 
Lya limit. The X-ray light curve shows an intense flaring activity 
at early times and it is dominated by the "late prompt" component 
at later times. The optical light curve is sampled only up to ^^2000 s 
rest frame (~ 9800 s observer frame) and is dominated by the "stan- 
dard afterglow" emission, as shown in Fig. [7] Since fluxes in only 
two optical bands are available it is not possible to infer the value 
of j4v°^* from the optical photometry. Curran et al. (2007) found a 
very soft observed spectrum, after correction for Galactic extinction 
and Lya absorption: the optical spectral index fi"^^ — 3.1 ± 0.4 
at 5000 s (observer frame). Assuming that the optical and X-ray 
emission are produced by the same mechanism they inferred two 
possible values of the host galaxy dust extinction, assuming either 
a single or a broken power-law joining the optical and the X-ray 
data. 

We extracted 2 SEDs at 4500 and 6500 s (observer frame; see 
Fig. |7). In both SEDs the X-ray data are extracted from the time 
interval reported in Tab. 1 . In the two component modelling the op- 
tical and X-ray emission would be due to different components, to 
account for their different temporal behaviours. Thus A^'^^ can- 
not be inferred from the optical to X-rays SED. This argument, 
together with the fact that only 2 optical bands have been sampled, 
does not allow a direct dust absorption estimate, therefore we do 
not have constraints also on the optical spectral index. If the opti- 
cal emission is produced by a standard afterglow mechanism, we 
can choose as an example a value of /3o ~ 0.5. If we assumed a 
(3o ~ 0.5 then A^'^^ « 0.65. Since only two optical bands are 
available, the uncertainties are very large but since this value is 
similar to the mean Small Magellanic Cloud-like -Eb-v = 0.27 
obtained by Curran et al. (2007) in the broken power-law case the 
latter value is used in the correction applied to the plotted SEDs. 
The large error on /9x,i, and the paucity of photometric data do not 
allow to draw any firm conclusion on this burst. 
The large uncertainties on the optical to X-ray SED at 4500 and 
6500 s (observer frame), and the choice of the same Eb-v = 0.27 
used by Curran et al. (2007) make them consistent with an unique 
broken-power law but the presence of two components cannot be 
excluded. 

No break is observed in the X-ray light curve. In particular. 
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Figure 7. Top panel: X-ray (in grey) and optical (different symbols, as la- 
belled) light curves of GRB 060210 (in the rest frame time). Same notation 
as in Fig. \5\ The vertical line and stripes indicate the jet times expected if 
the burst followed the -Epcak vs. Ej "Ghirlanda relation" (Ghirlanda et al. 
2007) (see text). The vertical lines mark the time at which the SEDs are 
extracted. Mid (Bottom) panel: Optical to X-ray uF^ SED of GRB 060210 
at about 4500 (6500) s after trigger in the observer frame (916 (1300) s rest 
frame). The dashed and solid Hues show respectively the best fit (with the 
ABP model) to the X-ray spectrum (the spectral parameters are reported in 
Tab. 2) and the uncertainties on the slope of the low energy spectral index 

/3x 1- The dotted line shows the best lit (with the AP model) to the X-ray 
^ ' , (c) 0000 RAS, MNRAS WO, 000-000 
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one can estimate the expected jet break time if the GRB was to 
follow the so called "Ghiiianda relation" (Ghirlanda et al. 2007). 
However the lack of evidence for such a break is consistent with the 
light-curve modelling as the X-ray flux is indeed dominated by the 
"late prompt" emission at the time when the jet break is expected. 
No data are available at such a time in the optical band where such 
a break should have been detectable, due to the dominance of the 
"standard afterglow" component. 

In conclusion, the poor optical photometry of GRB 060210 
does not allow to obtain an estimate of both A^^^ and /3o. This 
fact, together with the quite large errors in the /3x,i estimate, does 
not allow to infer firm conclusions on this GRB. The diversity of 
the optical and X-ray temporal behaviour and the lack of jet break 
in the late time X-ray observations prompt us to model the light 
curves as due to different components. The optical to X-ray SEDs 
cannot give better constraints to the model since within errors are 
consistent with both having one or two separate components. 

4.4 GRB 060729 

The UVOT data in 6 filters are from Grupe et al. (2007) while the 
ROSTE R band photometry is from Rykoff et al. (2009). 

After a steep decay in X-rays for about 400 s, the optical and 
X-rays light curves track each other and are characterised by a long 
lasting (~ 50 ks) shallow decay phase. Following an achromatic 
break, a steeper flux decay phase is observed in X-rays up to about 
10^ s after the burst (Grupe et al. 2009). The light-curves are domi- 
nated by the "late prompt" component with no significant evolution 
of the optical to X-ray flux ratio as can be seen in Fig. [8] As a con- 
sequence, the optical to X-ray SED is not expected to evolve in 
time. 

The bottom panel of Fig.[8]shows that the optical flux at about 
1500 s is indeed consistent with the extrapolation to the optical 
band of the broken power law X-ray spectrum. The SED does not 
require any additional host galaxy dust absorption and the poorly 
constrained optical spectral index is consistent within errors with 
/3x,i- The quality of the optical-UV SED is not good enough to 
directly constrain the Av°°'. 

It would be possible to consider SEDs at later times based on 
the UVOT data (the R band photometry covers only the first XRT 
orbit). However since neither the optical to X-ray flux ratio nor the 
colour significantly evolve we present here only the most complete 
SED at 1500 s. In this case the plotted X-ray spectrum is extracted 
from the first two PC mode orbits excluding the first 150 s of the 
first orbit in order to avoid the contribution of the steep decay phase. 



4.5 GRB 061007 

The photometric data are from Mundell et al. (2007) (/, R, V and 
B bands) and Rykoff et al. (2009) (ROTSE R band). After a steeper 
flux decay lasting about 90 s (rest frame), the X-ray light curve 
declines following a single power-law for the whole observed time. 
The optical flux instead shows a fast rise ( by about 2 orders of 
magnitude) in the first 40 s followed by a simple power-law decay 
up to about 60 ks rest frame (see Fig.[9}. The first R band fluxes are 
simultaneous with the 7-ray prompt emission, and the rise between 
the first and second detection is faster than f, hardy explainable 
with any standard emission mechanism. 

In the two component modelling the optical light curve (after 
the end of the prompt phase, i.e. ~ 50s in the rest frame) is domi- 
nated by the "standard afterglow" emission. The single power-law 
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Figure 8. Top panel: X-ray (in grey) and optical (different symbols, as 
labelled) light curves of GRB 060729 in rest frame time. Same notation as 
in Fig.|5] The vertical line represents the time at which the SED is extracted 
Bottom panel: Optical to X-ray i/Fu SED at about 1500 s after trigger in 
the observer frame (970 s rest frame). The dashed and solid lines show 
respectively the best fit (with the ABP model) to the X-ray spectrum (the 
spectral parameters are reported in Tab. 2) and the uncertainties on the slope 
of the low energy spectral index /3x i- The dotted line shows the best fit 
(with the AP model) to the X-ray spectrum. 



X-ray decay phase after the end of the 7-ray detection would be 
dominated by the "late prompt" component. Thus the X-ray and 
optical fluxes would be dominated by two different components, 
despite the similarity of the light curves after ~100 s (rest frame), 
requiring a hard "late prompt" /3x,i in order for this emission not 
to significantly contribute to the observed optical flux. 

The optical fluxes have been corrected for a host galaxy dust 
extinction = 0.54 ± 0.30 (Kann et al. 2009). We consid- 

ered two SEDs at the times where all of the four photometric bands 
are simultaneously available: the first one at about 270 s (observed 
frame), immediately after the beginning of the simple power-law 
X-ray decay, and the second one at about 5.5 ks (observed frame). 
The X-ray data in the latter SED are extracted from the first two 
orbits in PC mode while the ones plotted in the first SED are ex- 
tracted from the time interval reported in Tab. 1 . In both cases the 
hard /3x,i found with the broken power law fitting (see Tab. 2) im- 
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Figure 9. Top panel: X-ray and optical liglit curves of GRB 061007 in rest 
frame time. Same notation as in Fig. [5] The vertical line and stripes indicate 
the jet times expected if the burst followed the i?peak vs. "Ghirlanda 
relation" (Ghirlanda et al. 2007) (see text). The vertical lines mark the times 
at which the SEDs are extracted. Middle and bottom panels: Optical to X- 
ray vFv SED at about 270 s (middle) and 5500 s (bottom) after trigger in 
the observer frame (corresponding to 120 s and 2.4 ks in the rest frame). The 
dashed and solid lines show respectively the best fit (with the ABP model) 
to the X-ray spectrum (the spectral parameters are reported in Tab. 2) and 
the uncertainties on the slope of the low energy spectral index /3x i. The 
dotted line shows the best fit (with the AP model) to the X-ray spectrum. 



plies a negligible contribution of the X-ray component in the op- 
tical band, supporting the proposed interpretation as can be seen 
from the middle and bottom panels of Fig. [9] 

The X-ray light curve does not show any slope variation in 
correspondence to the expected jet break time obtained in the as- 
sumption that the GRB follows the "Ghirlanda relation". Once 
again this is in agreement with the "late prompt" dominated na- 
ture of the X-ray flux. The jet break should instead be visible in 
the optical, but unfortunately there are no observations after 150 ks 
(observer frame) to confirm or rule out this prediction. 

The early time optical to X-ray SEDs of GRB 061007 has 
been analysed also by Schady et al. (2007b) and by Mundell et al. 
(2007). They extracted the SED around 600 and 300 s observer 
frame after trigger respectively and they found these SEDs to be 
well fitted by a single power-law. As can be seen in fig. |9l the opti- 
cal fluxes at 270 s (observer frame) are consistent with an extrapo- 
lation of a single power law fit of the X-ray spectrum. Their single 
power law X-ray fits give results consistent with the ones presented 
in tab. 13.11 In this paper we consider also the broad band SED of 
GRB 061007 at later times. The bottom panel of fig. [9] shows that 
after 5 ks the optical fluxes are no more consistent with an extrapo- 
lation of the X-ray data single power law fit and the two component 
scenario that we considered in the light curve modelling is in good 
agreement with both the early and late time SEDs. 

4.6 GRB 061126 

A very rich photometric sampling is available for GRB 061126 
(Perley et al. 2008, Gomboc et al. 2008). After a steeper decay 
the X-ray light curve follows a single power-law flux decay for 
the whole observed time. The IR-optical~UV light curve instead 
shows a more complex behaviour, as shown in Fig. 1101 

We modelled the power-law decay of the X-ray light curve 
as a "late prompt" component. The optical bands is accounted for 
by a sudden transition from a "standard afterglow" dominated early 
time to a "late prompt" dominated late time behaviour. If correct, 
this scenario would imply a spectral evolution from a two compo- 
nents to a single component optical-to-X-ray SED. We would also 
expect an evolution of /3o at the time of transition between the two 
components. 

We extracted two SEDs, plotted in the middle and bottom pan- 
els of Fig. [To] The first one (middle panel) corresponds to ~2000 
s (observer frame) and is obtained using 8 contemporaneous pho- 
tometric bands (f/, B, V, R, I, J, H, and Ks). The X-ray data 
are extracted from the time interval reported in Tab. 1. At this 
time the optical spectrum is well fitted by a single power law with 
Po = 0.94 ± 0.05, and no host galaxy dust absorption is required. 
This is in agreement with the findings by Perley et al. (2008). The 
optical spectrum slope is inconsistent with the X-ray spectrum. We 
examined a second SED at 45 ks (observer frame; bottom panel of 
Fig. 1101 ) when both the optical and X-ray light curves are domi- 
nated by the "late prompt" component. The plotted X-ray data are 
extracted from the time integrated spectrum of the second XRT ob- 
servation. At that time only 4 optical bands are available (Ic, Rc, V 
and B), and the spectrum is still well fitted by a single power-law 
without host galaxy dust absorption, but the spectral index is harder 
than at earlier times (i.e. j3o = 0.54 ± 0.1). The bottom panel of 
Fig. 1 101 reveals that not only the optical fluxes but also the optical 
slope are now consistent with the extrapolation from the X-rays 
(with slope Px,i)- 

As predicted by the light curve modelling, the optical spectral 
index evolves after the transition from standard afterglow to late 
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Figure 10. Top panel: X-ray (in grey) and optical (different symbols, as 
labelled) light curves of GRB 061126 in rest frame time. Same notation 
as in Fig. |5] The vertical line and stripes indicate the jet times expected if 
the burst followed the iJpcak vs. E-y "Ghirlanda relation" (Ghirlanda et al. 
2007) (see text). The vertical lines mark the times at which the SEDs are 
extracted. Middle and bottom panels: Optical to X-ray vFv SED at about 
2200 s (middle) and 45 ks (bottom) after trigger in the observer frame (cor- 
responding to 1000 s and 21 ks in the rest frame). The dashed and solid lines 
show respectively the best fit (with the ABP model) to the X-ray spectrum 
(the spectral parameters are reported in Tab. 2) and the uncertainties on the 

tral index /3x i ■ The dotted line shows the best 
''"O-OOO ' 
ray spectrum. 



prompt emission, witli tlie SED becoming consistent witli a single 
dominating component. 

We can contrast our interpretation witii tlie alternative one pro- 
posed by Gomboc et al. (2008), who suggested that the presence 
of some dust absorption at early times could account for the opti- 
cal spectrum being consistent with a broken power-law optical-to- 
X-ray SED. At later times the SED could be fitted with a similar 
broken power-law but without the need of any host galaxy dust ab- 
sorption. Thus in the Gomboc et al. (2008) scenario a change in 
dust absorption would be required to interpret the optical spectral 
and broad band SED evolution. 

Also for this burst we can evaluate the expected jet break time 
in the hypothesis that the GRB follows the "Ghirlanda relation". 
However, at the corresponding epoch both the optical and X-ray 
emission are dominated by the late prompt contribution and thus 
no jet break would be observable. This is indeed in agreement with 
the absence of a break in the observed light curves, although the 
observations end soon after the predicted jet break time. 



5 CONCLUSIONS 

We analysed the Swift XRT data of a sample of 33 long GRBs se- 
lected by G09 to have known redshift, published estimate of the 
host galaxy dust absorption and good XRT and optical follow up. 

If the XRT 0.3-10 keV spectra are modelled as a single power 
law, we confirm that the host frame 7Vh°^* column densities are 
rather large when compared to the values of the host galaxy dust 
absorption inferred from the optical analysis, according to "stan- 
dard" extinction laws (see also Schady et al. 2007, Stratta et al. 
2004). 

For the 15 brightest bursts we could model the X-ray data 
with a broken power law, and in 7 cases we find evidence of a 
spectral break (90% confidence level). In such cases the required 
N^'^^ is in turn smaller than for the single power-law fitting and is 
marginally consistent with the column estimated by the optical ex- 
tinction. However, in other 8 bright GRBs the X-ray spectrum does 
not show any break and some of them do require a large value of 
jyhost -pi^gi-gfoi-e the presence of an intrinsic curvature in the spec- 
trum cannot be considered as a general solution for the "excess" of 
N^^^ commonly found in GRB X-ray spectral analysis. 

In order to test the interpretation by G09 that the X-ray and 
optical light curve complex behaviour can be interpreted as due to 
the contributions to the emission by two different components, we 
combine the results of the light curve de-convolution with the X- 
ray and broad band spectral properties at different times. 

In particular we checked whether the presence of a break in the 
XRT spectra is consistent with what has been observed in the opti- 
cal bands by studying the time dependent optical to X-rays SEDs 
of the GRBs for which a spectral break was found. 

We found that 7 of the 8 events are consistent with the pres- 
ence of a break in the XRT spectra and the evolution of the broad 
band SEDs appears to support the predictions of the two compo- 
nents scenario, even in the presence of complex light curve be- 
haviours. In one case (i.e. GRB 060210), the quality of the data 
does not allow us to solve the ambiguity between the temporal and 
spectral analysis. 

Consistency is also found in relation with the (lack of) evi- 
dence for jet breaks in the light curves, whose break time is esti- 
mated by assuming that the Ghirlanda relation holds for all GRBs 
(Ghirlanda et al. 2007). Indeed light curves are observed to steepen 
in correspondence with the jet break time only when the light curve 
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is dominated by the "standard afterglow" emission while no break 
is detected if the other component is dominating. 

Further testing of the two-component modelling requires to 
extend the simultaneous multi-band optical follow-up at later 
times (i.e. several days after the trigger), when the typical ex- 
pected R band magnitudes are around 24-25. Such an intensive 
and long lasting multi-band follow up could allow us: i) to search 
for possible optical spectral index evolution when a different com- 
ponent becomes dominant (as in the case of GRB 061126); ii) to 
test for the presence (absence) of jet breaks in "standard afterglow" 
("late prompt") dominated optical light-curves in a larger sample 
of events, presumably shedding light also on the jet geometry and 
energetics. 
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